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RESEARCH SUMMARY 

Major stand^'level and forest-wide decisions facing forest 
managers interested in applying even-aged management are identified. 
Traditional, current, or potentially useful techniques for making 
these decisions are reviewed. Emphasis is placed upon how each 
method characterizes the stand or forest, the growth and yield 
information requirements of each method, the comxDutational burden 
associated with each method, and the flexibility of each method 
and its adaptability to the form and complexity of the solution. 
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INTRODUCTION 


Effective planning and conduct of a forest management program requires the forest 
managei to make a number of basic decisions; therefore a thorougli understanding of the 
majoi questions prefacing these decisions is critical in forest management planning. 
Once the specific questions or decisions have been defined and understood, the forest 
manager chooses the particular analytical tool or tools to help in decisionmaking. 
Because a great deal of management science research has been devoted to developing 
alternative analytical tools, deciding which is best for a particular need can be 
difficult. The purpose of this study, therefore, is to examine the basic questions 
acing tlie manager interested in practicing even-aged forest management and to ex- 
plain and evaluate the various analytical tools historically, currently, or potentially 
useful in answering the questions. Evaluation criteria will include: the way a parti- 
cular method characterizes the stand or forest; the growth and yield information 
requirements of the method; the computational burden associated \tfith the method; the 
flexibility of the method as to the form and complexity of the solution; and other 
pertinent considerations. A similar study has recently been reported for uneven-aged 
management (Mann and Bare 1979) , 


Management Questions 

The basic timber management question is simply: What treatment should be applied 
to each stand each year to best meet the objectives established for the forest? The 
term ''treatment’' includes inaction. 


This apparently simple question is actually very complex and difficult to answer. 
The manager is faced with a forest composed of perhaps hundreds of stands; a myriad 
of potential treatment types, intensities, and timings from which to choose for each 
stand; and also with an increasingly complex set of objectives and constraints for the 
forest. Even in the early days of forest management when forestry objectives were 
simpler and the set of possible treatments was relatively simple (that is, harvest cut 
or not), finding an answer to the basic question was difficult. 

The basic question historically has been divided into a number of subordinate 
questions perhaps with the belief that answering subordinate questions muld be easier. 
These subordinate questions can be separated into two categories: stand- (or aggregrated 
stand-) level questions and forest-level questions. 

For stands (or stand classes) , the manager interested in practicing even-aged 
management needs to know the optimal: 

1 . Species mix 

2. Planting density 

3. Thinning plan 

4. Fertilization plan and 

5. Rotation length. 

While it has not been done historically, 2, 3, and 5 from above could be combined into 
the question, ’’IVliat is the optimal diameter or size distribution schedule?” 


1 



The term '’stand classes” refers to the practice of aggregating stands of similar 
species mix, site quality, and age into classes to reduce the coiiipexity of the problem. 
An ’’optimal” solution to a specified problem is defined to be that solution which best, 
or most favorably, meets the specified objectives and constraints of management. As 
defined, a plan includes both the schedule and intensity of treatment. 


At the forest level, the forest manager interested in applying even^aged management 
is faced with determining, for the desired number of planning periods, the optimal; 


1. Schedule of stand (or stand class) treatments and/or 

2. Conversion strategy and conversion period length. 

The two questions are closely related because a conversion strategy is also a schedule 
of stand treatments. The second question differs from the first if a targeted forest 
structure exists. Then the problem becomes how to best convert the present forest 
structure to the targeted structure while meeting the objectives and constraints of 
management. A regulated forest is an example of a targeted forest structure. Selection 
among alternatives for meeting one or several targeted forest structures constitutes 
the conversion strategy decision. The body of methods used are known as harvest 
scheduling techniques. Solution of these problems provides the allowable cut for each 
planning period. 

The above questions can be regarded as basic questions as well as the set of prob- 
lems dealt with in this paper. The techniques discussed in this paper can also be 
applied to make such subsidiary decisions as; 


1. Selection of the most efficient logging techniques to achieve specific stand 
objectives (Sessions 1979) 

2. Selection of the optimal harvesting schedule and of stand protection treatments 
in the face of such damaging agents as forest pests (Brodie and Rose 1975) 

3. Selection of the optimal schedule over time for both stand treatment and 
reading access (Navon 1976). 


These applications depend upon 
harvest sclieduling questions. 


a thorough understanding of the basic stand treatment and 
They will be elaborated on further in this paper. 


TRADITIONAL SOLUTION METHODS 

which have been studied the longest have been rotation length and 
the scheduling of stand class treatments to determine the allowable cut. ^ 

Traditional Rotation Length Methods 

As stated by Davis (1966), numerous criteria have been proposed to determine 

tlon,“t »i:?ci°a,s;a;rirrir«nfvT'^ regenera. 

this’age will the “JySS H TT' , 

culmination occurs will vary depending unon ^ stand. The age at which 

stand and upon the measure of vSIume^ 

pressed in total stem cubic foot volume than it sooner when ex- 

(Davis 1966). ^ expressed in board foot volume 



Maximum forest rent is conceptually identical to culmination of NJAI , the only 
difference being that tlie latter maximizes average annual volume growth while the former 
maximizes average annual net value growth. For a given stand age, annual net value 
growth is the gross value of the timber minus growing costs all divided by the stand 
age (Davis 1966). As an economic criterion for determining rotation length, this method 
has been strongly criticized for ignoring the opportunity cost incurred by locking up 
capital in a forest investment. 


The final method, maximum soil (or land) expectation value, does include the 
opportunity cost of capital. The first step in computing rotation length using maximum 
soil expectation value is to determine, for a given rotation length and interest rate 
the present net worth of all costs and returns, incurred and/or received during the 
rotation. Soil expectation value of the specific rotation length is then computed as 
the present value from receiving an infinite series of these present net worth "payments. 
This process is repeated for a number of rotation lengths and the rotation with the 
largest soil expectation value is chosen. The rotation of maximum soil expectation 
value will provide the oiimer with the highest return on investment in the 'property at 
the specified discount rate (Davis 1966). 


Evaluation of Traditional Rotation Length Methods 

CHARACTERIZATION OF THE STAND 


Traditional methods require the stand to be characterized by tliose factors 
to obtain access to the particular yield tables used in the metliod's analysis, 
actors usually include site quality, stand age, and species. 


needed 

These 


GROWTH AND YIELD INFORMATION 


in "i^thods use yield information, the assumptions incorporated 

the yield information will subsequently be used in determining the rotation length. 

unmanaged stands are used, it is assumed that 
the stand will be fully stocked throughout its life span and that no mortality capture 
y inning will occur. On the other hand, if the yield tables incorporate specified 
thinning plans then it is assumed that the thinnings will occur on the specified 

specified intensity. While all methods do determine an optimal 
sn?MHnn specific criteria, it should not be assumed that the 

^ represents the optimal treatment program for the stand. Rather, the 

solution IS optimal for the given treatment program incorporated in the yield tables. 

COMPUTATIONAL BURDEN 


^ All traditional methods can be done by hand. The product size 
simplest method to apply, while the maximum soil expectation metliod 
information and is the most difficult to compute. 


criterion is the 
requires the most 


FLEXIBILITY OF SOLUTION SPECIFICATIONS 


Each traditional method has been designed to meet a specific objective. ' 
objectye of the product size criterion is to produce trees of a specified, fi- 
The objective of the maximum mean annual increment criterion is to produce the' 
amount of wood volume possible, while maximum forest rent has the objective of 
the greatest average annual net value of the stand for a single rotation. Fin; 
objective of the maximum soil expectation criterion is to produce the greatest 
net value of the stand for an infinite number of rotations. 



Traditional Methods for Scheduling Stand Treatments 
to Determine Allowable Cut 

of .tllln f^stion that has historically received much attention is the scheduling 

find solutions J^hirm^l allowable cut. Three methods have been used to 

lutions to this pioblem-volume control, area control, and area-volume check. 

aoDlv^^ volume control requires the least amount of information to 

a information is required and so allocation of ^h^ 

pitted cut to actual stands must be accomplished after the fart TUf^. n t -r 
control is to provide a constant volu,e of cutover a spe^Uied oerlS 
have been developed to compute the allowable cut in this fashinn <;n ' 4 ^ 

IrSrErS 

number of acres cut llch yL is siiDlvMirto?ar?^"°t length for the type, the 
In years. The usual rule for alloratinn .i forest area divided by rotation length 

( 0 ” «anfn:;s,"L\1”h'e llZf 

C«hic"is'',nofrSt‘ej'tjrcase? “applSLVof^""* P«'“'«ivity classes 

traditional approach o tSI LSl™ ha bLn difficult. One 

type classes il, equivalent prSiScJir^v un?,; to express the different productivity and 
fractional number of acres upprlpri tr. ^ j * units are usually expressed as the 
harvest. iSereforf d'riihod f in regulated final 

control (Davis 1966). Electronic comT)utpr 5 control instead of true area 

(or productivity unit) control to theso simplified the application of areas 

and Chappelle 1967) is an exanple of an area ?«pi“'^proJS!'“ ““ 

can Piovidrf!omsta!ralloMbl'rM,“u“;oMra«as' for'”' control .ethods 

harvested in each planning period. Regulated forest^Jr ^'nits) being 

encountered, which leads to nrohlemc ofests are, however, not usually 

control results in an erratic number of acres^f^'^'^ unregulated forest, volume 

produces an erratic volumrJut The IZTf "! harvested while area control 
control formulas will be only a crude annroximaHon o”® the volume 

can be sustained during conversion to a^lnTsti f n .p optimal harvest that 

third methodology of scheduling stand clLs^t^eat^e^tf to 1"^^® problems, a 

evolved, the area-volume check method. ^eatments to determine allowable cut has 


which try to combine feSurerfrom^both^th°'^^^? ^ techniques, all of 

Because of this diversity it is diffiLu ■"'® "’"^hods. 

methods do, however, require detiilpd info ?®’^®ralize about the procedures. All 
stand Clashes, and asLTthaf th^^lver™^ or 

methods also require a great deal of tedious cLculatiorto length. Most 

figure. The last feature has made the introdStioJ of ?hrA?f . ® “ allowable cut 

area-volume check problems a blessing to many foJSte^s! ®^®''^’'onic computer to solve 
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As an example of area-volume check, the following is a brief description of a 

S^euriShhl '' of the USDA Forest Service 

I I -f- -4-1 * first Step in this method is to obtain an initial estimate of 

stLrfJlaS dat^So^r ‘J*’ formulas such as Hantlik's formula. The 

flr^; , rr'i ArF AiTai'Sed by cutting priority, which is usually "oldest stands 
f; t , IT required to harvest each stand class is then calculated using th? 

? ^ ^ ^ expected yield from each stand class at final harvest is cal 

vest rhrfnr^^r first examines whether the actual time needed to completely har- 

IcrL^e h^riesLr^th conversion period and the second examinL whether the 

^I'erSe aUowaMe cut these checks are unacceptable. 

ChaLelL fT!. "" systematic fashion and the process is repeated, 

perfori the opeiations^'"'’ computer program ARVOL to actually 


Evaluation of Traditional Stand Treatment Scheduling Methods 

CHARACTERIZATION OF THE FOREST ° 


Volume control through formula does not need to characterize the forest because it 

chSactLizrjhe fo?esfh^®i°^ treatments. The other two traditional methods 

HAflt ^ K ^ forest by lumping stands into rather broad classes that are usually 
defined by age and perhaps site quality. usuaiiy 

GROWTH AND YIELD INFORMATION 

rvAA fheir yields are fixed for each application of the method 

Except for repeated solutions, no choice of how a stand class can be treated is provided- 
for example runs with and without thinning might be compared. In addition, thinning 
u s can on y be expressed as increases in final harvests or as priority harvests fthat 
, thinning harvests must be taken before regeneration harvests) . 

COMPUl’ATIONAL BURDEN 


All of the traditional methods can be computed by hand, 
ormula is the easiest to compute, while the area-volume chec 
clitricult and tedious technique to apply. 


Volume control through 
method is tlie most 


FLEXIBILITY OF THE SOLUTION SPECIFICATION 


None of the traditional methods 
scheduling. 


can handl,e economic criteria for stand treatment 


OTHER CONSIDERATIONS 


Both area control and area-volume check assume that the conversion period equals 
the rotation length and that a targeted forest structure exists. In addition ^h^ 

uS'‘SiS“rc';:t°oide:r rr tL 

harJeS priority! determine the 


5 



OPERATIONS RESEARCH METHODS 

The introduction of both the electronic computer and operations research techniques 
grciitly expanded the manager's ability to answer questions pertaining to even-aged stands 
and forests* Sofne of the first applications of the computer were to mechanize the 
computations of numerical values previously, and tediously, done by hand, such as the 
aforementioned allowable cut methodologies. Two broad areas of operations research have 
seen application to forestry problems, simulation and optimization. 

Ihe distinction between simulation techniques and optimization techniques in even- 
aged stand and forest-level analysis is not distinct because of confusion in usage. 
Confusion arises from a body of intermediate methods lumped together under the term 
binary search techniques," We will distinguish between pure simulation, pure opti- 
mization, and binary search techniques. For the purpose of categorization, we will 

however, follow the usage of Johnson and Tedder^ and lump the binary search techniciues 
under simulation. ' 


/ pure simulation approach requires that the analyst specify the entire treatment 
p an tor each stand. The simulation tool used provides a detailed report of the impacts 
these decisions. Sensitivity to changes in these decisions is provided by repeated 
thrie^ro/^n^Tf treatment plans. The alternative select^ wmild be 

constraints ^ tried, judged against some implied objectives and set of 

constiaints. Supeiior alternatives could easily exist, but may be missed. 

pna I 'T® ^Pq^i^ation approach, the analyst supplies only a specific obi active 

to finrthfbLt^leverof the optimization technique then procLds 

tinn ^ ^ ^ objective criterion that can be met under the snecifica 

the .Of course;- 

sJecIficatTonf lijl'opSmiSiVappr f t^"" the p?obLi'unJer"a?teJmate 

mathematical programing techniwes SclSa°"d ^^en-aged management have used 
integer programing. These applications will ’be^disJuLe^’^'"’'’ nonlinear, goal, and 

approach. ^Advantage can^b^Lken^of knwledee^af ^ implicit in the pure simulation 
problems to efficiently solvrsoL about the structure of forest management 

binary search techniquL'ni anLer lienfLa"? ’1710 

under certain objective criteria and constrainr^ ^ ^ optimization tecliniqucs 

solutions under others. Currently binary searrli^t different and lower valued 
in handling complex sets of regional inveMor^dita ^ ^^f®^ provide greater flexibility 
forest tyi)es with different allowable cut assumptions?'^°^'^^"^ multiple ownerships and 


■ I * I w I I 

Simulation Methods for Answering Stand-Level Questions 

1 - 01 3 ta„d- 

the user to modify the thinning plan and the^rotatior? ^ manner which allows 

possible consequences of the specified 

B nr program. Each simulation run produces 

^Johnson, K, Norman, and Phil in i t u 
1" <=« calculation. Por.'^Sci, 
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information concerning stand development, yields from thinning and harvest cuts, and 
present net worth based on discounted costs and revenue figures. This information is 
reviewed by the user (the forest manager) to see if the particular solution looks "best” 
when compared to other selected management programs. 


addition to examining thinning plans and rotation lengths, Hamilton and Christie 
.. ® with an opportunity to examine initial planting densities, while 

Hoyer (1975) provides the means of examining the effects of fertilization at one 
intensity. 


Ev^ludtion of Simuldtion Methods for Answering 
Stand-Level Questions 

CHARACTERIZATION OF THE STAND 


The particular data items used to characterize the stand depend upon the type of 
stand development model used by the simulation method. These data items can range from 
site index, stand age, and stand basal area and/or total number of trees to site index 
stand age, and the size and competitive position of individual trees within the stand. ^ 

GROWTH AND YIELD INFORMATION 

Stand simulation methods can be applied to any kind of existing stand development 
model. With the use of more current stand development models, it is often possible to 
predict yields from a wide array of treatment plans. 

COMPUTATIONAL BURDEN 


Most stand simulators require a computer to operate. 

FLEXIBILITY OF THE SOLUTION SPECIFICATION 

While the stand simulation methods can be used to assess either physical or economic 
objectives, the methods cannot assure that the examined management program represents 
the global optimal solution to any specified set of complex objectives, constraints 
and cost-revenue relationships. To determine a global or near global optimal solution 
would require a great number of simulations and careful analysis and can be accomplished 
more efficiently by optimization for many stand development models. 

Simulation Methods for Answering Forest-wide Questions 

Simulation has also been used extensively in answering forest-level questions, 
even before the introduction of the electronic computer. The area-volume check method 
IS an example of a precomputer binary search simulation method. The 
course, increased the manager's capability to utilize simulation. It 
the size and complexity of problems that can be addressed with ■ 
shortened the time and tedium associated with obtaining a solu^ 
examples of where binary search simulation has been employed f 
managerial questions are the programs: SORAC (Chappell' ' ~ 

(Sassaman and others 1972), and TREES (Tedder and oth- 
application of pure simulation is program TEVAP2 (M^ 

Myers 1974) . 

Program SORAC can compute an allowable cut us 
volume check. The program allows the manager to a 
period in the planning horizon. If area control i 
parameters, then the solution is identical to that 
control parameters that change between periods wil 
differ from AREA, When the area-volume check meth 
horizon, the results from SORAC also differ from t 
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Progi’ain ARVOL computes an allowable cut that can be sustained for the conversion 
or rotation period. Using this method, however, can cause the allowable cut to change 
drastically after the conversion period or first rotation. Also, a manager rarely 
computes an allowable cut that will be followed for an entire conversion or rotation 
period. Instead, he will apply the cut during the current planning period and recompute 
a new allowable cut figure at the staz't of the next planning period. Using this continual 
updating feature, drastic changes in allowable cut from one planning period to the next 
can be avoided. SORAC basically tries to simulate this actual planning and updating 
process. While the allowable cut for the current planning period is the same as that 
provided by ARVOL, the advantage of SORAC over ARVOL is that SORAC provides the manager 
with a more realistic picture of the long-term consequences of his actions. 


Like ARVOL and AREA, SORAC characterizes the forest by stand classes that are 
usually based on stand age. Finally, either extensively or intensively managed yields 
can be specified for each stand class, with intensive yields Implying a predetermined 
treatment schedule. 


In program SIMAC, the forest is also characterized by stand age classes, with a 
user-specified average site index for the whole forest. SIMAC, however, allows much 
more flexibility in specifiying a treatment scheme for testing than SORAC. SIMAC also 
ditfep through its allowable cut procedures. First, no rotation or conversion period 
length IS specified. Instead, the length of the planning horizon is specified fup to 
40 decades) along with the youngest age in which a stand class can be regeneration 

lu allowable cut, the simulation can harvest a stand class at any 

age do™ to the minimum age (that is, final harvest age is not predetermined). 

difference is that an age-related cutting priority can be specified 
instead of assumed. A third difference is that a fixed allowable cut can be specified 

cu c^ Z bound on aHowirie 

cut can be specified and the SIMAC program will compute, if possible, the maximum 
wUhin ihes^bS^f planning horizon and still be 

e.bodyJ„V?l;f ‘r"’ i' “ techniques 

approaches (ECHO is discussed under onf J control, SIMAC, SORAC, or RCliO type 

1, Harvest priority; 

2. Variable silvicultural intensity over time and portions of the inventory^ 

sllvISuJral »ea»ra?! '"’’"sl'ed. ownership, species, site, and prescribed 

presently available hS''resultS”in*tSsrtechniquerb ''l'” eptii»izatlon packages 

projections (Beuter and others 1976). ^ chosen for detailed regional 
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The final example of an allowable cut simulator is program TEVAP2, Besides being 
a pure ratlier than binary search method, this program has a number of features which 
distinguish it from SORAC, SIMAC, and TREES, First, TEVAP2 deals with stands rather 
than stand classes. The stands must be even-aged with one or two stories, but they can 
be uniquely characterized by their type; site index; and the average diameter, height, 
number of trees, age, and dwarf mistletoe rating of the overstory and understory. For 
this reason, average or ’’normal” stands are not assumed and treatments can be better 
tailored to the actual stand. It should be noted that Myer's (1968, 1973) stand 
simulators can easily be used to derive the treatment schedule to be used in this forest - 
level analysis. 

The second distinguishing feature of TEVAP2 is its planning horizon. In TEVAP2, 
the planning horizon is equal to the planning period and, for many users, tlie planning 
period is one decade, TEVAP2, therefore, is a short-range planning tool when compared 
to SORAC, SIMAC, and TREES, 

Allowable cut determination methodology is the final distinguishing feature of 
TEVAP2, Four allowable cut levels are computed for the planning period. The first 
allowable cut figure is calculated under the assumption that all specified treatments 
are applied to optimally stocked stands existing in a forest with a perfectly balanced 
age structure, The second allowable cut is calculated under the assumption that all 
specified treatments are applied to the existing stands and forest structure. The 
third cut is computed by applying area regulation to the existing forest structure. 

Tlie final allowable cut value is computed by volume control using a modified version of 
the Austrian formula. 


The allowable cut figures are not applied to the existing forest structure for 
the conversion or rotation length to determine if they can be sustained, Feather, these 
cuts serve as guides and the forest manager can redefine treatment schedules and inten- 
sities if he is not pleased with the results. 


Evaluation of Simulation Methods for 

Questions 


CHARACTERIZATION OF THE FOREST 


Answering Forest-wide 


Some forest -wide simulation methods (for example, SORAC and SIMAC) characterize 
the forest by combining stands into classes using stand age, site quality, and/or 
species. Other forest-wide simulation methods (for example, TEVAP2 and TREES) charac 
terize and maintain each stand of the forest using detailed stand information. 

GROWTH AND YIELD INFORMATION 


Some forest-wide simulation models, such as SORAC, SIMAC, and TREES, use tabular 
yield information which can be derived from any yield source. Other forest-'"^ 
simulation models, such as TEVAP2, have stand development models built into 
stand development models usually have been of the whole stand type. Unlike 

traditional forest-wide decision methods, all of the aforementioned fore„_ 

simulation methods, except SORAC, allow thinning yields to be credited to total harvest 
when they occur instead of at the final regeneration harvest cut. 

COMPUTATIONAL BURDEN 

All of the forest-wide simulation methods used as examples require a computer to 
operate. IVlien compared to forest-wide optimization techniques, simulation methods are 
often faster to run on the computer and therefore cheaper (Johnson and Tedder^). 
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FLEXIBILITY OF THE SOLUTION SPECIFICATION 

Simulation methods SORAC, SIMAC, and TEVAP2 are used exclusively for seeking the 
objective of physical wood production. Program TREES can be used to meet either 
physical or economic objectives. As with the stand- level simulators, the treatment 
schedules and intensities are fixed for the run while cutting priority may be fixed by 
the simulation method. To determine a truly global, optimal solution would require 
repeated trials with different schedules, intensities, and priorities. If the potential 
number of treatment tyi)es is large, including such types as planting intensity, thinning, 
fertilization, harvesting method, and species mix, finding a global solution may be 
extremely difficult. 


Another limiting feature of the forest -wide simulation methods is the difficulty of 
introducing constraints into the analysis that might be applied to the liarvest flow. 
Simulation models, however, can more readily identify feasible solutions than their 
optimization counterparts (Johnson and Tedder^). 


OPTIMIZATION METHODS AT THE STAND LEVEL 

Optimization methods for solving stand-level questions have long been used in 
orestry. Early methods of determining rotation length, such as maximum MAI and soil 
expectation were all optimization methods. These early methods were distinguished by 
the numerous assumptions made in order to solve them and by the amount of tedious 
calculations required. First applications of the electronic computer were aimed at 
so ving tie attei problem. More recently, mathematical programing techniques have 
made it easier to solve more complex problems and, therefore, reduce the need for 
numerous assumptions. A good example of this evolutionary process has been the develop- 
qu^t ^ optimization tools for jointly answering thinning plan and rotation length 


Optimal Thinning Plan and Rotation Length 

of thinnine^nlan”flnH^^^^"«- presented an early solution to the joint optimization 

aL JaJS^fl analvsisTol^S^v equations 

net profi? LraiZ rut l °P levels which would maximi ze 

net protit per annum. Then given an initial stocking level, the ontimal stockina 1ovr.i 

'''®>' determined the volL to be removed by tMmSg in ‘ 

tional load. First, the growth models werne^y'^mplistif ^ Second° 

was assumed to be either constant ov^r ^ * npiistic. Second, the stunipage price 

age. Pinally, costs that “"rie3 .irrJLHl' “f ^taml 

method (thatch, stcpeg* :r:ssu“li llTZT.ltZl 

could be used to obtain"'the’same*^so\utions^L^r?^^^ demonstrated how dynamic programing 
programing is a discrete m^tSaticaf^ ra^^^ ""'f" 

through a network of points. The numbL of a ^ describe a problem 

point also defines the dimensions of the network needed to define each 

each point by its volume stocking, with 1 000 hoard -F J^idon and Akin defined 

by its stand age, with 5-year in^^rvaK LS? intervals between points, and 

two state descriptors and was two dimensional*" pomts. Their network, therefore, had 
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The objective of dynamic programing is to find, within defined limits, the best or 
optimal path through the network. In the Aiiiidon and Akin (1968) problem, therefore, 
the objective was to find the optimal stocking level at each age class. Vhe limit s’ of 
the network are defined by physical considerations such as maximum growth of the stand 
and the impossibility of negative volumes. Solution to the problem of finding an 
optimal path can be accomplished using either the forward or the backvvard recursive 
method , ^ 

Because the forward recursive method is easier to understand, it will be explained 
fiist using tlie Ainidon and Akin problem as an example. An initial stand age and stocking 
level aie specified which defines the starting point in the network. A growth function 
is then used to advance the stand to one of the points in the next time period. If 
nothing were done the stand would remain at the stocking level of the point to which it 
had grown. Thinning, however, can reduce stocking which shifts the stand to other 
lower points in that time period. For the second time period, there is only one path 
the stand could have followed to get to each point because of the single starting point 
in the first time period. Because the second period has defined a number of points that 
the stand can grow from, there are a number of different paths in which a stand could 
arrive at a stocking point in the third time period, Sfliat is wanted, though, is the 
optimal path to that point. Therefore, all possible paths from one age to the next are 
examined and only information for the optimal path is stored for each possible stocking 

level at time period three. This process is continued for all stand age classes in the 
network. 

The values that must be stored at every point in the network, excluding the starting 
point, are the value of the objective function for the best path to that point and the 
coordinates of the previous point on the optimal path. With the forward recursion 
method, therefore, it is possible to examine all possible points at a given age and, 
using the objective function value, determine the optimal stocking level for that age. 
Given that point, tlie optimal path (and therefore stocking level for each age) can be 
determined from the stored coordinates. 

Inherent in this method of determining the optimal path is the "principle of 
optimality." Translated into the problem of determining the optimal thinning plan, 
this principl e states that, once the optimal thinning plan has been determined to a 
specified stand age and structure, the optimal plan for the next older stand age depends 
only on the older stand's age-structure combinations not yet analyzed. Thus, the various 
possible stand structures for a given age need to be analyzed only once, which greatly 
reduces the number of calculations necessary. 

The backward recursion method ultimately provides the same answers as the forward 
recursion method, but the values stored at each point are different. The backward 
method starts at all possible ending points and determines the best nath to tho startino 
point. Stored at each point, therefore, is the optimal path to 
value of the objective function. 

The advantage of the backward recursion ir'" 
from the optimal path due to insects, fire, or 
already computed will show which new path is o 
age. The disadvantage (and therefore the adv" 
is that a new run must be made for each potei 


^It is our belief that dynamic programii 
in the decisionmaking process. We also belie 
difficult tool to understand for many in our 
proportionally more attention to it. 
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Aniidon and Akin (1968) used the backward recursive method. They also greatly 
restricted the problem by allowing the stand to increase or decrease only 1,000 board 
feet in one growth period (this is, between adjacent age periods). This meant that a 
given point could only be reached from two adjacent points. While this reduced the 
complexity of the problem and made solution easier, it also greatly restricted possible 
thinning plans. 


Because the backward recursive method solves for only one rotation length in a 
given run, Amidon and Akin (1968) had to make a number of runs to evaluate the effect 
^ notation length. Like Chappelle and Nelson (1964), tliey used the maximum 
soil expectation criterion for determining the optimal rotation length. 

As pointed out by Amidon and Akin (1968), the advantage of dynamic programing over 
the marginal analysis approach of Chappelle::and Nelson (1964) is its computational 
flexibility. Once the problem has been structured for dynamic programing, it is easy 
to change pricing and cost assumptions. Also, the results of one run can be used to 
evaluate other "suboptimal" solutions that might become viable if stand coiiditions 
change. A final computational advantage is that stand development models and dynamic 
functions need not be as simple as those used by Chappelle and Nelson (1964) and Amidon 
and Akin (1968). 


After examining the approaches of Chappelle and Nelson (1964) and Amidon and Akin 
(1968), Schreuder (1971) made two critical comments on their methodology. First, he 
contended the four had ignored cost of land when determining optimal stocking levels. 

This omission, he said, can affect the determinination of economic stocking levels. 
Second, he felt that their method would not find the jointly optimal thinning plan and 
rotation length because their method "...does not consider the possible interdependencies 
of the period." Therefore, their method would ignore thinnings with negative value 
for the period but which, over the rotation period, could more than pay for themselves. 
Their method would also ignore timber selling costs, which vary with volume removed. 

Both of these factors could lead to suboptimal paths. 


Schreuder (1971) proposed that thinning schedule, intensity, and rotation length 
could be easily optimized jointly if the final harvest cut were defined as an extreme 
inning, e formulated the joint optimization problem in a continuous, control- 
theoretic form in which revenue, costs, rate of timber cut, stand yield, and growth 
were all fimctions of time. This formulation would determine the optimal thinning and 
^antity of volume removed in continuous, instead of discrete, values of time and volume, 
^ce he had the formulation, however, he discovered that an explicit solution could not 
be found except for trival problems. Schreuder (1971) then turned to the numerical 
solution method of dynamic programing to find a discrete solution to his formulation. 
Although he formulated the solution in a backward recursion, dynamic programing frame- 
work, he did not show actual examples of the process, 

f programing to solve the joint optimization problem 

has recently been reported by Brodie and others (1978) and by Brodie and Kao (1979). 

In addition, Brodie and others (1978) demonstrated that, contrary to Schreuder 's (1971) 
findings, the approaches of Chappelle and Nelson (1964) and Amidon and Akin (1968) do 
indeed simultaneously solve for both optimal stocking level and rotation length. 

The purpose of the Brodie and others (1978) study was to assess the "impacts of 
regeneration cost, initial stocking level, site, quality premiums, and variable logging 
costs on thinning and rotation...." Three features distinguish this dynamic programing 

affaS! aJ Ji “ «omUty estimator. More realistic estimates 

evaluated^ ^ ^ thinning, through mortality capture, therefore, can be made and 
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Second, the dynamic network of Brodie and others (1978) allows much greater accuracy 
in defining thinning intensity than did Aniidon and Akin^s (1968) network. From any 
specified node, the latter’s network allowed the stocking level to either increase or 
decrease 1,000 board feet. Brodie and others (1978) allow stocking to decrease by 100 
cubic foot increments to zero and to increase by 100 cubic foot increments to maximum 
potential stockings for the age class. 

Finally, Brodie and others (1978) chose to use the forward recursion method instead 
of the backward as advocated by Amidon and Akin (1968) and Schreuder (1971) . For 
thinning analysis, they considered the forward method a more flexible tool than tlie 
backward method. 

One shortcoming that Brodie and others (1978) identified with their approach was 
that the stand model did not provide diameter growth acceleration usually expected from 
thinning. If quality premiums are important, this shortcoming can lead to suboptimal 
solutions. 

The most recent work in the application of dynamic programing to stand problems, 
by Brodie and Kao (1979), eliminates this shortcoming. They accomplished this by 
using a much more complex stand model that incorporates growth acceleration in the 
quadratic mean stand diameter. After a careful examination of the stand model, Brodie 
and Kao (1979) found that the model could be initialized for a specified site index 
if the three values of stand age, basal area, and number of trees were known. These 
variables, therefore, formed the three state descriptors for the dynamic programing 
solution. The sizes of the intervals between nodes are 10 years, 4 square feet of basal 
area, and 15 trees, respectively. ' 


Analysis results using the three-descriptor framework are the optimal number of 
trees and basal area to maintain in each time interval. To obtain these values, Brodie 
and Kao (1979) used the forward recursion solution method. Because quadratic mean 
stand diameter can be computed from number of trees and stand basal area, stumpage 
prices and logging costs could be more realistically introduced into the anaylsis as 
functions of quadratic mean stand diameter. 

Other Possible Methods of Determining Optimal Thinning Plan 

and Rotation Length 

The dynamic programing approaches reported to date have all ignored or assumed a 
diameter distribution structure. Diameter distribution data is important j however, in 
the planning of milling facilities and as an aid in applying specified treatments to 
field conditions. Therefore, a procedure that can determine the optimal number of 
trees for each time period by diameter classes may have greater usefulness than present 
procedures . 

Conceptually, dynamic programing can also be used to determine optima'’ ' 

class distributions over time. The state descriptors ’ ’ 

(T^) in one of the Dth diameter classes at one of th" 
the number of dimensions for the netvvork space is D+ 

network is A(T^) . From these formulas, it is easy t 
fast becomes impossible in practical applications. E 
to determine, within 10 trees per diameter class, the 
tain in each 2-inch diameter class at the start of ea 
that the maximum number of trees per diameter class w 
120 years, and the maximum diameter class size was 4C 
classes would be 26 (0, 1-10, 11-20, ..., 241-250), t 
120), and the number of diameter classes 20 (2, 4, .. 
number of dimensions for the network space would be c 

of nodes in the network would be 12x26^^= 3.12x10^^. 
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Discussion in this section has been confined to deterministic analysis of stand- 
level questions. Stochastic stand-level decision analysis using dynamic programing is 
presented in several papers (Hool 1966, Lembersky and Johnson 1975; Leinbersky 1976). 

Another potentially promising approach has been described by Adams and Ek (1974). 

Part of their work dealt with the problem of converting an existing uneven-aged diameter 
distribution to a targeted structure over a number of time periods in which cutting 
could occur. This problem is analagous to the even-aged situation of starting with 
a specified stand, growing, and thinning it over a number of time periods until a targeted 
diameter distribution of zero trees in each diameter class is reached in the last time 
period, Adams and Ek (1974) used a nonlinear programing approach which necessitated a 
relatively simple stand model that explicitly predicted the net change in number of 
trees in a diameter class. Their method was also limited as to the number of time 
periods over which the method could be applied. 

It is obvious that methods for determining what diameter distributions are optimal 
to carry over time are in their infancy. Alternatives to dynamic or nonlinear program- 
ing, such as application of optimal control theory or decomposition theory, have been 
suggested (Adams 1974 and Adams and Ek 1976) , but the development and tesing of practical 
applications is limited. 


Optimal Fertilization Plan and Species Mix 

Fertilization and species mix are two additional factors which can interact and 
affect optimal thinning plans and rotation lengths. Consequently, they should all be 
jointly optimized to reflect their interdependencies. It has already been shown how 
an optimal thinning plan, and rotation length can be jointly determined for a stand. 
The pioblem now is the incorporation of fertilization and species mix into tlie process 


The optimization of qming and intensity for fertilization has not been addressed 
in the literature. This is explained by the scarcity of stand models with built-in 
fertilization response functions. Once appropriate stand models have been developed, 
the joint optimization of thinning schedule and intensity, rotation length, fertilization 
schedule, and intensity can be done using dynamic programing. Incorporation of fertili- 
za ion intensity results only in a mulitiplicative increase in computations and not an 
exponentiation of network size, 3 provided that entry-time intervals are longer than the 
effect of fertilization. 


Species mix is another area where little work has been reported. As with fertili- 

snecies “’"c Jen"rmn?t simulators that can handle multiple 
species. Given a multiple species stand model, optimization could be done by use of 

?rti;rp?oHr'"'’ descriptor would be required for each specLs 


Evaluation of Optimization Methods for 
Answering Stand-Level Questions 

CHARACTERIZATION OF THE STAND 


As with the stand-level simulation methods, all stand-levi 
characterize the stand by those features used to initialize th( 
development model type used by the optimization methods. 


Optimization methods 
particular stand 


practices--S™ercia!^and prccoMLSal'’tSninr’fertir silvicultural 

effort. PH.p. diss. »ep.^f for. ita»^;s:;™irrani:^'"sfr 
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GROWTH AND YIBLD INFORMATION 


Stand “level decision using dynamics programing can be applied to any type o£ stand 
development model, if it is acceptable to the manager to prespecify thinning rules 
(such as, cut smallest or largest trees first;) for priority of removals (Martin 1978), 

If this condition is not acceptable, then application of dynamic programing is currently 
limited to the whole-stand/diameter^free type of stand development model. 

Analyses using nonlinear programing, and perhaps control theoretic techniques, are 
potentially applicable to the whole-stand/diameter^class type of stand development model. 

COMPUTATIONAL BURDEN 

In application of dynamic programing, size of the problem can cause difficulties 
and pose limitations. For example, if a stand simulator existed that was initialized 
by number of trees and basal area for any of a number of species and if A is the number 
of age classes, T is the number of tree classes, B is the number of basal area classes, 
and S is the number of species, then the dimensional size of the dynamics programing 

c 

network space would be 2S+1 and the number of nodes would be A*(B«T) . Thus, problems 
of this type become large very fast and, if a diameter class distribution type of 
simulator were used, the size of the problem is increased many fold (that is, using 
these and the previously defined symbols, size of the network space would be SD+1 and 

S • D 

number of nodes would be A*T * ). 

c ^ 

Another limitation of optimization methods is the size of their computer programs. 
This problem is particularly aggravated if the optimization method is interfaced with 
some of the complex stand development model t>qDes. 

FLEXIBILITY OF THE SOLUTION SPECIFICATION 

Both the dynamic programing and the nonlinear programing methods can determine the 
optimal treatment plan to follow to meet both physical and economic management objectives. 

OTHER CONSIDERATIONS 

As a precaution, the ’’principle of optimality” that is inherent in dynamic pro- 
graming is a restriction that must be dealt with carefully. It has been previously 
shown how a model that characterizes stand development by site index, age, basal area, 
and number of trees per acre can be used to determine optimal thinning plans by the 
dynamic programing technique. Because the time variant variables (age, basal area, 
and number of trees) were all state descriptors, the ’'principle of optimality” can be 
met. If additional variables are used in the model, then one of three actions can be 
taken to still meet the principle. First, the additional variables can be added to the 
analyses as additional state descriptors. This is the most obvious and theoretically 
best solution, but, as demonstrated, it comes at the potentially high c ' " 

requirement of larger computer capacity. 

The second method of handling additional variables is to make assui 
the relationship between the additional variable (s) and the state descr: 
was the method used by Martin (1978). Finally, for some variables, it : 
avoid the problem by modifying the existing state descriptors. Suppose 
that additional variables (such as time since last cutting or time sinc) 
have been used to explain variation not explained by the state descript< 
residual basal area and/or number of trees), If the effect of these va: 
zero after a period of time, then these variables can be handled, withoi 
additional state descriptors, by defining the planning period length to 
long as or longer than the period it takes for the effect to disappear. 
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OPTIMIZATION METHODS AT THE FOREST LEVEL 

Ovei' the past two decades, a considerable amount of effort has been directed at 
developing techniques for answering forest-wide management questions. An excellent 
evaluation, analysis, and summary of most of these techniques has been provided by 
.lohnson and Schcurman (1977) . They defined two basic mathematical programing formula- 
tions, which they called Model I and Model II, and two pricing options. The two basic 
fonnulations differ by their treatment of the stand classes. For the sake of simplicity, 
.Johnson and Scheurman (1977) assume that the forest is composed of one species and one 
site quality aiid that the resulting forest structure can be adequately portrayed by age 


With this forest structiu'e, Model I defines a management unit by its initial age 
c ass and these units are maintained throughout the analyses. Therefore, only one set 
0 activities consisting of a number of alternative treatment plans must be defined for 
each management unit. A treatment plan is defined by the timing of harvest and thinning 

returns expected across the entire planning horizon. The optimal 

i-rf'iMn ' by the number of acres in a management unit assigned to each 

tieatracnt plan in order to meet management objectives and constraints. 

<1odtl II also defines starting management units by their initial age classes New 
management units, however, are defined if harvesting of the initia age^classero^curs 
K. IS. If harvesting occurs in several different age classes durinf a period thel 

j.os5il,lc timings of the fitst ImrvLt ’(,.(" 0 ^“^ ar(he'OTd’'of“(!'^ ““ 

,n:ur’;(:t“s„r;;::n"“ra“rdof^ 

timings of the second han'est cut or of the roti harvesting all possible 

if the units arc never reharvosted The rcmT the end of the planning horizon 

type is egnal to the ^rter o^ircujrJSt"™'^* of activities of this 

horizon to a single stand class. planning 


and Scl;i‘.™r,;‘(l™5)^';‘vfsl’o™''i,'‘rtie "rae't’ivm^s ' 

Area constraints guarantee that the area harvested^! Model I formulation, 
an age class in any period. NonnegaJIvStJ JoSstr^?c"° within 

and areas from a class are not negative. ^Other coJst^ainrr"^®^ volumes 

control constraints, can also be incornorated in rw' harvest volume 

objectives. The number of harvest control constr^^nr!’'^^^^^! managerial 

tion The preceding emphasis on turnumter orac;i 1 
modcl to formulate a solution is because of the' concJJr?ha? 

do not exceed computer limitations. ^ ^ resulting formulations 
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Linear Programing 

As mentioned earlier^ each formulation for solving the optimum scheduling of stands 
in the forest can have two pricing options if the objective involves an economic 
criterion. The first, and historically the most widely used, option is founded on the 
principle that for a given stand class, treatment plan, and planning period, the unit 
price received is constant regardless of the volume cut. If this assumption is 
acceptable, the operations research technique of linear programing can be applied to 
solve both the Model I and the Model II formulations. The works of Loucks (1964), Navon 
(1971), and Ware and Clutter (1971) represent the application of linear programing to 
solving Model I formulations. In their review, Johnson and Scheurman (1977) demonstrated 
how the simulation techniques of the area-volume check method, SORAC (Chappelle 1966) 
and SIMAC (Sassaman and others 1972) could all be reformulated into different forms of 
Model II. They also showed how Nautiyal and Pearse's (1967) formulation for solving 
the optimal even-aged conversion problem could also be formulated as another form of 
Model II, 


Mathematical Programing Techniques for Quadratic 

Objective Functions 

The second pricing option assumes that for a given stand class, treatment plan, 
and planning period, the unit price received for an acre cut is a function of the volume 
cut, and the form is one that declines linearly as the acreage cut increases. If this 
assumption is acceptable and reasonable, then, theoretically, the operations research 
technique quadratic programing can be applied to solve both basic formulations. A 
limited and specialized algorithm form has also been developed that incorporates the 
quadratic objective function. The first example of this algorithm form was Walker's 
(1976) ECHO (Economic Harvest Optimization Model), which Johnson and Scheurman (1977) 
reported to be a special formulation of Model II. A revised form of this algorithm 
is incorporated in TREES (Tedder and others 1980). Hrubes and Navon (1976) have 
approximated the ECHO solution using separable linear programing. 

Goal Programing 

Both Model I and Model II pertain to management programs that are concerned with 
single objectives such as maximization of volume production or of present net worth. 

For programs that are concerned with multiple objectives, the historical approach has 
been to pick one of them for optimization and treat the rest as constraints. In the 
area of multiple objective optimization techniques, goal programing has received some 
attention in forest related problems. Kao and Brodie (1979) have looked into the use 
of goal programing for scheduling the treatment of all stands within the forest. They 
examined the problem of maximizing present net worth (objective 1), of maintaining an 
even-flow of harvest (objective 2), and of producing a regulated forest structure at 
the end of the planning horizon (objective 3). In another study. Field and others*' 
have demonstrated the trade off of maximizing volume while maintaining nondec lining 
evenflow on a southern forest. 


**Field, Richard C., Peter E. Dress, and James C. Fortson 
and goal programing procedures for timber harvest scheduling. 
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Goal programing is an extension o£ linear programing. The objective of goal 
programing is to minimize the weighted difference between the targeted value for each 
objective and the feasible value of each objective. The weights of priorities assigned 
to each objective are specified by the analyst and can be changed between runs to 
evaluate their influence upon the solution. Using a simple example, Kao and Brodie 
(1979] demonstrated that goal programing could produce reasonable, compromise solutions 
to problems with infeasible linear pi'ograming solutions due to conflicting objectives. 

Integer Programing 

Both linear and goal programing can result in the acreage of a stand class being 
assigned to a number of alternative treatment plans. If it is desired that a stand 
class be assigned only one treatment plan, the appropriate optimization tool is either 
integer or integer-goal programing. To date, the computational limitations of these 
tools have precluded their widespread application to practical forestry problems (Bare 
and Norman 1969). Recent developments in optimization technology (Glover and others 
1978), however, may provide for wider use of integer programing in the future. 


Evaluation of Optimization Methods for 
Answering Forest-wide Questions 

CHARACTERIZATION OF THE FOREST 

In the forest -wide optimization methods, forests are usually characterized by very 
simple structures based on gross stand class aggregates. Aggregating stands in this 
fashion could cause problems in accurately predicting treatment plan responses because 
stand development models are usually nonlinear in form, 


The disadvantage of having to aggregate stands has risen out of the need to minimize 
problem size so that solutions are computationally feasible, that is, the fewer stand 
classes, the smaller the problem. There is, however, another feasibility problem which 
may conflict with the objective of reducing problem size by defining few stand classes. 
IVhen constrained optimization is used, it is often found that a large number of alternate 
treatment plans must be defined for each stand class in order to find an operationally 
feasible solution. 


As an alternative to this historical approach, it seems equally possible to define 
a large number of more homogenous stands (or if need be, stand classes) and a few 
alternate treatment plans for each. These alternate treatment plans could be developed 
by determining the optimal treatment plan for several different rotation lengths. The 
choice of rotation lengths could be made through an intelligent examination of existing 
forest structure problems and by determining where resulting future yields may be most 
lacking. This approach would produce more defensible yield predictions, provide a 
link between stand level and forest-wide optimization, and would tie treatments to the 
basic management unit, the stand. Testing this approach for feasibility seems justified. 

Problems associated with this disadvantage may be further minimized (or even 
eliminated) by the recent development of a new method for solving linear programing 
problems (Kolata 1979). This new technique is reportedly capable of solving much 
larger problems in shorter time periods than the currently used simplex method. 

GROWTH AND YIELD INFORMATION 

The yield information for all of the forest -wide optimization methods is tabular 
in form.^ Therefore, any source of yield information can be utilized. As with the 
forest -wide simulation methods, thinning yields can be credited to the total harvest 
when they occur instead of at the final regeneration harvest cut. 
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COMPUTATIONAL BURDEN 


Forest -wide optimization methods usually take longer computer times to operate 
than their simulation counterparts (Johnson and Tedder^)* Recent advances (Kolata 1979), 
however, may minimize or eliminate tl\is disadvantage, 

FLEXIBILITY OF THE SOLUTION SPECIFICATION 

All forest-wide optimization methods can evaluate both physical and economic 
objectives. Some of the advantages of optimization methods include: 

1 . Alternative treatment plans for each stand class can be evaluated in one run 
and an optimal plan determined, 

2, Cutting priorities do not need to be specified nor assumed. 

3, Complex constraint specifications can be incorporated into the analysis. 

4. Solutions can usually be generalized because first order optimization conditions 
are stated or implied in the analysis method. 

Finding feasible solutions, however, is often more difficult with optimization methods 
than with simulation methods (Johnson and Tedder^), 

OTHER CONSIDERATIONS 

In many of the optimization methods, fractionization of stands can occur by 
assigning more than one treatment plan to the stand. The use of integer programing is 
a possible future solution to this problem, 

SIMULTANEOUS, OPTIMAL SOLUTION OF THE 
STAND- AND FOREST-LEVEL QUESTIONS 

If a perfectly regulated forest structure exists or if harvest-level constraints 
are not important, then the joint optimization of stand and forest-level questions can 
be computed easily by determining each stand's optimal treatment plan and aggregating 
across the whole forest. Unfortunately, these conditions seldom exist and the manager 
is left with the problem of trying to determine suboptimal stand treatment plans that 
will meet forest -wide objectives and constraints. 

Work on this joint optimization problem has only recently received some attention, 
Williams (1976) and Nazareth (1973) independently demonstrated an approach for joint 
optimization using an extremely simple forest structure. They found that it is the 
size of the joint optimization problem and not its structure that leads to solution 
difficulties. Many of the linear-programing harvest scheduling studies discussed above 
provide for optimization over several stand- level treatments. The dynamic programing 
solutions discussed above indicate a large, finite number of stand-level treatments. 

A single, large linear programing formulation, incorporating all of these stand- level 
treatments as activities, would be unmanageable in size. The problem is solvable, 
however, by decomposition--a process by which large linear programing problems can be 
partitioned into smaller problems whose solutions can be aggregated to provide the 
solution to the orginal large problem. 
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For the joint stand-level and harvest scheduling problem, the process is described 
as follows: 

1. Solve the dynamic programing stand-level treatment network to find the optimal 
stand -level treatment sequence and the optimal sequence of treatments for every 
other potential stand condition, A large number of potential stand treatment 
activities are thus defined. 

2, Fnter the optimal sequence of stand-level treatments as activities in a con- 
strained linear programing harvest scheduling problem and solve for the 
optimum harvest schedule. This solution is usually not the optimum harvest 
schedule over all possible stand-level treatments as there are probably non- 
optimal stand-level treatments that enhance the constrained objective function. 
For example, heavier than optimal thinning or early> final harvest can provide 
volume in periods where available volumes from the optimal stand solutions are 
constraining, 

3. A decomposition is performed whereby new activities are selected using infor- 
mation from the first solution as a starting point, A new optimum harvest 
schedule is formed using these activities, 

4, Successive decompositions are performed with the process terminating after a 
prespecified number of decompositions or when the objective function enhance- 
ment between decompositions is trivial, 


Because this joint optimization technique combines both dynamic and linear pro- 
graming, it also exhibits their strengths and weaknesses, Tliis teclmique lias been 
used to jointly determine the optimal thinning plans, rotation lengths, and the harvest- 
ing schedule for a forest of simple structure and a single species (Nazareth 1973; 
Williams 1976). While theoretically feasible, incorporation of optimal fertilization 
plans, species mixes, and diameter distribution for more complex forest structures has 

The recent advance in linear programing technology (Kolata 
1979} may, however, eliminate the need for decomposition to solve this or even more 
complex problems. In any case, the linear programing portion of the method still would 
cause the fractionalization of stands (classes). 


SUMMARY 

The assignment of the best (or optimal) treatment plan to every stand in the forest 

y, this basic decision was made by partitioning it into several separate decisions 

d^^sj^r'thim'nrsr^an^'"' subordinate decisions have included the'opuLl plaJung 
density, thinning plan (and resulting diameter distributions), rotation length 

scE^aiC'rf'stand 

to ha^dirsiLly o?iSntlv th^ decisionmaking tools have been developed 

might be possible using nonlinear progrLinr JhP VJ®®® particular decisions 

from existing tools appears alquarS tSf’tS^ P°^^ible 

optimal diameter distributions,^fertilLation nlan^ Stand-level decisions concerning 
not received adequate attention. Growth LdelinrilmiJat^^®''''®^ however, have 

factors holding back development in area! ® be the major 
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At the forest level, a number of tools have also been developed to help in making 
stand treatment scheduling decisions. Additional research on the usage of integer 
programing to schedule entire stands and on the development of multiple objective 
optimization techniques would provide even greater analytical flexibility to the 
decisionmaker. 

Finally, a technique for the simultaneous consideration of both stand--level and 
forest“Wide decisions has recently been developed and reported (Nazareth 1973; Williams 
1976), Expansion of this technology to include optimal diameter distributions, fertili- 
zation plans and species mixes for multiple objectives should ultimately lead to the 
solution of the basic managerial decision. 
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GLOSSARY 

“f tl'« can bn removed in a planning period if the 

optimal schedule of stand treatments is fully applied to the forLt. 

Allocation problem - The task of distributing limited resources to alternative uses. 

programing formuiatlons to limit 

r ' real area constraints is necessary to limit area harvested 

from an age class to an amount less than or equal to the arL in thf agrclas^ 

Area control - A classical harvest scheduling technique that achieves a regulated 

0? the ;«S^d' -=1’ tli^p^L * i orient 

?ong-^: SuS-JurduHn'”?,' V har.osting areas of constant 

depen5ro^rt\fSi‘!\rc'":rd"i?:"S2;io:! 

““ thrarif°coM:of::So7 -i"® 

Area-volume-check^-fAirepetitlve classical harvest scheduling technique that attempts 

forest An inHf 1 1 volume, area, and time period of harvesting over the 

volume estimate of harvest is applied to the forest arerand 
resulting areas harvSted. Length of Um^for 
estimate is if found unsatisfactory, the initial 

a« JudJers^lsSto^r’'''^ ““ 

to the youngest stand ti,» o! solves the optimal treatment plan 

if the stand should devfaie SL t7Jptl.af oatV^e^t'^'r^?" 

£^rhour?ovi“ 

next iteration trial until a^^Lre'Sru^^r^S^^tl^Slst^!'’” 
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accumulated or depleted to reach the fully regu^ed SndijLr'"'"'' 
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Decomposition (Dmitzig Wolfe decomposition) - A process by which Icirge linear programing 
problems can be disaggregated to provide tiie solution to the original large problem. 

Deterministic analysis ^ Analysis conducted under the assumption that variables are 
single valued (that is, not governed by a probability distribution). 

Diameter distribution scliedule - A scliedule for tlie number of trees by diameter class 
for each planning period. 

Dynamic programing (DP) - a technique for solving mathematical programing problems using 
discrete numerical computation rather than simultaneous solution of continuous 
functions. Used in forestry for simultaneous solution of thinning plan and 
rotation length, and for other allocation problems. 

liCHO (Economic //arvest optimization) - An algorithm (Walker 1976) that uses a binary 
searcl'i technique for maximizing the present net worth of harvest scliedule under 
the assumption of a downward sloping demand curve. The solution approximates a 
quadratic programing solution with specialized restrictions. 

Feasible solutions - Any solution that satisfies all objectives and constraints in a 
problem specification. The optimal solution is the maximum or minimum valued 
feasible solution. 

Fertilization plan - A program of fertilization that incorporates both the scheduling 
and the intensity of fertilization. 

Forward recursion - A dynamic programing solution strategy that starts with the youngest 
stand age and determines for each possible rotation length over the planning 
horizon the optimal treatment plan for each rotation length. The advantage of 
the forward recursion method is that optimal treatment plans for all rotation 
lengths are determined in one run. In addition, its method of handling the growth 
of stands makes it easier to adapt stand development models to a dynamic programing 
forinulat ion. 


Goal programing (GP) - A technique for solving mathematical programing problems with 
multiple objectives and goals by minimizing the deviations from the goals. Two 
uses in forestry have been to determine compromise solutions for both multiple 
uses or harvest scheduling objectives. 

Harvest scheduling technique - A method for determining when individual stands or stand 
classes should be thinned or final harvested in order to meet the forest owner’s 
objectives. 


Integer programing (IP) - A technique for solving mathematical programing problems where 
the variables are constrained to take integer values. Used in forestry for harvest 
scheduling problems where complete stands must be assigned to a particular treat- 
ment plan use. 

Linear programing (LP) - A technique for solving mathematical programing problems in 
which the objective and constraint variables are linear, Used in forestry for 
solving resource allocation and harvest scheduling problems. 

LP activities - A set of treatment plans that can be applied to a specified management 
unit . 
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LP constraints - Specified relationships that guarantee that a problem solution meets 
both the conditions desired (or imposed) by the manager and conditions necessary 
to meet the problem specification. (See Area constraints, Nonnegativity constraints, 
and Volume constraints for specific examples.) 

Management unit - A group of stands usually defined by their initial stand age. 

Marginal analysis - Analysis based on the economic optimization principle that marginal 
revenue should equal marginal cost, where marginal revenue is defined as tlie 
contribution to revenue of an additional unit produced and marginal cost is defined 
as the contribution to cost of an additional unit of production. So long as 
additional production contributes more to revenue than cost, production should be 
expanded. 


i fathemat ical progiaming ^ A body of optimization techniques for solving problems with 
unctionally specified objectives and constraints. The maximum or minimum value 
of the objective is found by a defined process in a finite number of steps. (See 
anear piogiaming. Goal programing. Integer programing, Nonlinear programing. 
Quadratic programing, and Dynamic programing for specific examples.) 


Maximum torest rent rotation - The rotation for which average annual revenue minus 

is greatest. As an economic criterion for determining rotation length, this 

TocSl. ul criticized for ignoring the opportunity cost incSr;ecrbJ 

locking up capital in forest investment. ^ 


Haxin,™ mi C«can annual Incre.er,t) Rotation - The rotation for wl.lcl, volume divided 

rim. rotation provides the greatest physical yield in the long- 


«a.,l.„. prosont - The rotation for .hid, discounted revenue einus 

soil e.xpectation rotation. rotation is longer than the maximum 

"'^'Tsc\tt:f:osrirgj::SreJ^^ -venue minus 

rotations. This cruSiofvawfi^ n I*'® basis of an infinite series of 
use for given species, prices costs'^ and value in even-aged forest 

than an evaluation ba^ed ora’sJngl^ro^atJSn 

"^"iL;tttro1^S:1:ltLTm^-^^^^^^ Which the 

and few eonstrnint, .hen compared with d,e 

units as they are harvested. Model II formulatimif^h^^^^#'^ management 

more constraints than the comparable Lri foCLS. "^"" activities and 

'lr^aphf^L?,';SrL^^^”E™ as two or more dimensional 

graphical space characterized by: volume and aee developing in a 

age, and basal area. The numbed of va^LS^ ?!; volume, 

the dimensions of the network. Dynamic nrna-r^L ^ clescriptor variables) define 
discrete intervals within the network. ^ analysis is undertaken using 
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Nonlinear programing (NLP) - A body of techniques for solving mathematical programing 
problems in which the objective and constraint variables are nonlinear. Used in 
forestry to determine uneven-aged diameter distribution schedules, conversion 
schedules for uneven-aged stands, and thinning plan and rotation age for even- 
aged forests. 

Nonnegativity constraints - Constraints imposed in mathematical programing formulations 
to limit variables to zero or positive values because negative volumes or acres 
harvested are infeasible and have no meaning in the usual problem specification. 

Normal yield tables - Tables that present, by site index, and stand age, the yield of 
^'fully stocked,” undisturbed natural stands. 

Optimal solution - A solution which best, or most favorably, meets the specified 
objectives and constraints of management. 

Planning horizon - The total number of planning periods in an analysis. 

Planning period - The shortest time period (usually in years) of concern in an analysis. 
One-, five-, and ten-year planning periods are common. 

Principle of optimality - The fundamental principle on which dynamic programing is 

based. The principle states that once the optimal path through the dynamic pro- 
graming network has been determined from the initial state to a specific state, 
the optimal path through the remaining stages of the network depends only on the 
portion of the network not yet analyzed, The principle of optimality allows a 
great reduction in the amount of necessary calculations. 

Pure simulation - A technique for examining the consequences of applying a set of 

management decisions and/or rules to a stand or forest. Alternative management 
decisions and/or rules can be examined through repeated application of the 
simulation process. 

Quadratic programing - A special case of nonlinear programing in which the objective 
function is a quadratic form. 

Regulated forest (even-aged) - Strictly defined, a forest with a stand age class distri- 
bution such that the area in each age class from age one to rotation age is equal 
to the total area of the forest divided by the rotation age. If this forest is 
harvested on an annual basis at the rate of area/rotation age, the yield and age- 
class structure will remain constant over time. The term in practice is applied 
to forest age-class distributions that are approximately uniformly balanced so 
as to approximate the above equilibrium conditions. 

Rotation length - The period of time elapsing between the time a stand is regenerated 

and the time it is final harvested. Applicable to the even-aged management system. 

SIMAC (simulated Jntensively Managed Allowable Cut) - An algorithm (Sassaman and others 
1972) for computing the allowable cut over the planning horizon using binary search 
techniques. Age related cutting priorities, fixed allowable cuts (for the first 
ten decades), and an upper and lower bound on the allowable cut can all be 
specified by the user. 

Single-tree/distance-dependent model - A stand development model that characterizes the 

stand by a plot composed of individual tree measurements, including tree coordinates. 
Using the tree coordinates, an individual tree’s competitive status within the 
plot is computed by examining the nearness and size of surrounding trees, Models 
of this type have been developed for even-aged and uneven-aged, pure and mixed 
species stands. 
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Single-tree/distance-independent model - A stand development model that characterizes 

the stand by a representative sample of individual tree measurements. Because tree 
coordinates are not used, a tree^s competitive status within the stand is charac- 
terized by comparing the tree's size to all other trees in the stand or sample. 
Models of this type have been developed for even- and uneven-aged, pure and mixed 
species stands. 

SOflAC (5hort Run Allowable Cut) - A binary search algorithm (Chappelle and Sassaman 
1968) for computing the allowable cut using area control, area-volume check, or 
both. The area control option provides greater flexibility than algorithm AREA, 
while the area-volume check method has incorporated into it a dynamic look-ahead 
feature which was not in ARVOL, The dynamic look-ahead feature simulates the 
continual updating process normally associated with timber management planning, 

Stand class - An aggregation of stands of similar species mix, site quality, and age 
into a common class. 


Stand development model - A mathematical model (often coded for computer use) tliat 
predicts the development of the tree component of a stand over time. Ideally, 
a stand development model should incorporate components for: gross growth, non- 
catastrophic mortality, regeneration, catastrophic losses, the introduction of 
managerial decisions, and the portrayal of the prediction result. (See Whole- 
stand/diameter-free, IVhole-stand/diameter-function, IVhole-stand/diameter-class, 

Single-tree/distance-independent models for 
definitions of the various types of stand development models.) 


State descriptors - A set of variables that uniquely characterize a condition (state) 

t iat a system may be in. It is used in dynamic programing analysis where even-aged 
stands are characterized by two-descriptors (age and volume or age and basal area) 
oi three descriptors (age, basal area, and number of trees) . 

Stochastic analysis - Analysis conducted under the assumption that variables can take 
distribSJf probability (that is. governed by a probability 


TEVAP2 (Timber Evaluation And Planning 2] - A pure simulation algorithm (Myers 1974-) 
for examining alternative allowable cut strategies over a single planning period 

incorporates both the schedultM and the 
TreatmentiPlan^-tA of^»reat.ent(s) that Incorporates both the scheduline and 

aSri'd Srf ooonC’ “d'- 

volume control, SIMAC SORAC or ECHO ^he alternatives of area control, 

other algorithiix; tbpcc n u ^XP® approaches. When compared to these 

.athe^lcal progra.l„e tarnations to 
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Volume control - A classical formula based harvest scheduling technique that uses 

forest volume growth, inventory or age class data in various combinations to deter- 
mine conversion period allowable harvest. Volume harvested under this method will 
be constant for the planning period, but area harvested will vary depending on the 
initial age-class distribution. Formula results are gross approximations which 
were usually adjusted using area-volume -check techniques. 

Wiole-stand/diameter-class model - A stand development model that provides average 
stand information by characterizing the stand as the number of trees in defined 
diameter classes (that is, dividing the stand diameter distribution into classes 
and predicting the number of trees in each class) . Models of this type can be 
developed for even- or uneven-aged, natural or managed stands of pure or mixed 
species. They are initialized by site index, number of trees by diameter classes, 
and, perhaps, by additional stand variables such as stand age. 

liniole-stand/dlaraeter-free model - A stand development model that provides average stand 
information, but does not provide diameter distribution data. Models of this type 
are usually for pure, managed or natural even-aged stands, and are initialized by 
site index and stand values such as stand age, average or quadratic diameter, and 
the basal area or total number of trees in the stand. 

Whole-stand/diameter-function model - A stand development model that provides average 
stand information through the use of smooth functions, such as tlie Weibull or beta 
distribution functions, to characterize the diameter distribution. Models of this 
type have usually been limited to pure unmanaged even-aged stands and are initialized 
by site index and stand age. 
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